Gravity compensation in complex plasmas by application of a temperature gradient 
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Micron sized particles are suspended or even lifted up in a gas by thermophoresis. This allows 
the study of many processes occurring in strongly coupled complex plasmas at the kinetic level in a 
relatively stress- free environment. First results are presented. The technique is also of interest for 
technological applications. 



Dusty or complex plasmas are low temperature, low 
pressure plasmas, e.g. glow discharge plasmas, contain- 
ing microparticles in addition to ions, electrons and neu- 
tral gas. They have a large variety of important applica- 
tions in fundamental and applied physics (see e.g. 0). 
For example, the microparticles can, under certain con- 
ditions, form a crystal |g, ^, ^ which allows the study of 
the liquid-solid phase transition at the kinetic level 
Due to the large mass of the microparticles (typically 
10^^ times the mass of plasma ions) gravity is an impor- 
tant force. The normal procedure for complex plasma 
experiments under gravity conditions is to suspend the 
microparticles (which carry a negative charge, Q) in a 
strong electrostatic field, E, i.e. mg = QE + ^^Fi, 
where the Fi are other forces such as ion drag. To low- 
est order, E varies linearly with position in a sheath of 
extent and it is easy to calculate that systematic body 
force variations from one lattice plane to the next exceed 
interparticle forces for the interesting particle size range 
above a few /^m. The systems are thus subject to a con- 
siderable amount of stress and there is abundant stored 
energy, which affects especially nonlinear processes such 
as phase transitions, interface dynamics and configura- 
tional or velocity space instabilities. It is for this reason 
that experiments under microgravity were proposed and 
are now being conducted 

It is the purpose of this paper to demonstrate that for 
monodisperse particles thermophoresis can be employed 
to counteract gravity, allowing complex plasma experi- 
ments to be performed in the "stress free" environment 
of the central nearly field-free plasma of a RF discharge, 
and to present first measurements. Results show that 
large stable systems can be produced in the charge neu- 
tral section of the plasma which exhibit the " central void" 
generally observed under 0-g condition. New is the ob- 
served saturation of thermophoretic levitation in the case 
of a large particle number combined with low working 
pressure. Jellum et al. were the first to apply a tem- 
perature gradient and to identify the thermophoretic ef- 
fect in a dusty plasma. In order to compensate gravity, 
however, more prerequisites have to be fulfilled: a homo- 
geneous temperature gradient by design, very low energy 
input to the gas, and the use of monodisperse particles. 

Using elementary kinetic arguments, the ther- 
mophoretic force can be estimated similar to the heat 



conductivity g. The latter follows from the energy flux 
in the presence of a temperature gradient. Instead of the 
energy transfer, we consider the momentum transfer per 
unit time and area in a plasma. Then the thermophoretic 
force on a microparticle with radius is given as 
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where v = {SkT/nm,)^^'^ is the average thermal velocity 
of the gas atoms with mass to at a gas temperature T 
and k is the Boltzmann constant. A is the coefficient of 
heat conductivity and dT / dx the temperature gradient. 

As we will see, this formula provides very good results 
if measured values or the best theoretical calculations are 
used for A. The following result for mono atomic gases 
has been given in the literature [|| 



A 



2.4 — , 

TO 



(2) 



where c = 3fc/2 is the specific heat per atom and rj the 
shear viscosity of the gas. For the viscosity, we use a 
result derived from a precise solution of the transport 
problem in the case of hard spheres B 
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with the gas kinetic cross section a for atomic scatter- 
ing, which can be taken from the literature jl^, Q . Us- 
ing these cross sections, values for the viscosity of noble 
gases are obtained (see Table I) in good agreement with 
measured ones. For example, the measured viscosity of 
Argon at T = 300 K is = 2.21 • 10"^ Nsm'^ compared 
to 2.16- 10^^ Nsm~^ following from (^ together with the 
cross section a — 4.2-10"^^ vo? given in Table I. Combin- 
ing (|l|), (H), and (HI) results in the thermophoretic force 
on a microparticle 
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Like viscosity and heat conduction the thermophoretic 
force is pressure independent, as can been seen from (Q) 
containing no pressure dependent parameters. Eq. (|4|) is 
valid only for spherical particles, mono atomic gases and 



2 



low pressure where the mean free path is much larger 
than the particle radius. 

We now compare our result to the literature. First, 
Waldmann jl^ derived the thermophoretic force using 
the Enskog-Chapman method for solving the Boltzmann 
equation. He obtained the same expression as (|^), where, 
however, the coefficient is given by 32/15 instead of 8/3, 
i.e., a thermophoretic force reduced by about 20% com- 
pared to our result. 

In order to compare our result with the one by Havnes 
et al. @, we use the relation nX = l/(\/2cr) g in (|), 
where n is the number density of the gas and A the mean 
free path of the atoms, leading to 

F = -4.67 nkXrl —. (5) 
^ ax 

The result obtained by Havnes et al. follows from re- 
placing the factor 4.67 by 8 in (|^), i.e., they found a 
thermophoretic force which is larger than ours by almost 
a factor of two. Balabanov et al. |Q have used a formula 
for the thermophoretic force in a complex DC plasma 
which agrees with (|) assuming the ideal gas relation for 
the pressure and replacing the coefficient 4.67 by 4. 

TABLE I: Gaskinetic atom-atom cross-sections a, mean free 
path A(p) at 50 Pa, viscosity and heat conductivity at 293 K, 
interaction radii R = 2r — \foJi^ from two sources. The 
thermophoretic force F is compiled for the actual particle 
radius — 1.69 fim and temperature gradient of 1170 K/m. 
For a given temperature difference, F is strongest for Helium. 
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Measurements were done in a completely symmet- 
rical RF-excited plasma chamber with a volume of 400 
cm^ (Fig- 0)- The gas is heated from below and cooled 
from above by Peltier elements. The RF amplitude (with 
frequency 13.56 MHz) is applied to the electrodes with 
180° phase difference. 

For calibration of the thermophoretic effect the dis- 
charge was operated at only 20 mW. About 10'* parti- 
cles were injected, much less than shown in e.g. FigJ^. 
The temperature difference between the metal plates 
((4) in Fig.|l|) was increased until two narrow particle 
clouds, above and below the central void (Fig.^), of equal 
size were visible. For Argon gas the equilibrium was 
reached for temperatures of 54.6 C at the lower elec- 
trode and 19.5 C at the upper electrode. The temper- 
ature difference corresponds to a gradient of 1170 Km-^. 



As monodisperse particles still have a residual size and 
weight spectrum, particles observed above the void are 
different from those below in that they are about 5% less 
heavy. The 50/50 partition hence calibrates the ther- 
mophoretic force by the average particle weight. For the 
measured temperature gradient of 1170 Km~^ (^ pre- 
dicts a thermophoretic force of 3.66 • 10""'^^ N. Our parti- 
cles (monodisperse melamine resin, diameter 3.4±0.1/im, 
density 1510 kg/m^) have a weight of 3.0 ± 0.3 • lO-^^ 
Thus, formula (^ is confirmed within the experimental 
uncertainties. We repeated the calibration with Neon 
and found equilibrium at a temperature difference of 
17.7 K. Thus thermophoresis in Neon is twice as strong as 
in Argon as expected from the gas kinetic cross-section in 
Table |. We also verified the 50/50 partition as described 
above for a pressure range of 8 to 62 Pa and found no 
pressure dependence in thermophoresis as predicted by 
theory. 




FIG. 1: Cross-section through the chamber: (1) stainless steel 
electrode, (2) dust dispenser, (3) Macor insulator, (4) metal 
plate at chassis potential, (5) window (welded Vycor body), 
(6) Peltier element, (7) heat sink, (8-10) fields of view (ther- 
mophoresis on: dashed line, off: dotted dashed) . The distance 
between the electrodes is 30 mm. 

For comparison and further description of the appara- 
tus we start with a condition where the Peltier elements 
and thermophoresis are not active (Fig. H). The particle 
cloud settles down above the lower electrode due to grav- 
ity. One observes that the particles form a lattice, the so- 
called plasma crystal. The plasma glow was filtered away 
by an interference filter peaked on the laser wavelength 
of 686 nm. Measurements taken without filter indicate 
the highest particle density congruent with the brightest 
glow. There is no convection after the particles have ar- 
ranged themselves into the ordered structure. The laser 
cut is focused to «50/im full width, much narrower than 
the lattice spacing. The stringlike structures, caused by 
the wake potential of the positive ions streaming to the 
lower electrode, disappear in patches, corresponding to 
local distortions no longer illuminated. The three di- 
mensional structure of such a system was analyzed by 
Zuzic et al. [|l^ . The RF power measured at the genera- 
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tor is 170 mW, the peak to peak amplitude measured by 
a voltage probe at the electrode is 45V. Both electrodes 
are shunted with 3.3 kfi to chassis. Losses in the unbal- 
anced/balanced transformer, the matching circuits and 
in particular the shunt resistors dictate that only 10% 
(17 mW) of the RF power go to the discharge. 




Electrode 



FIG. 2: Side view of 3.4/^m microparticles suspended in 
a plasma (complex plasma), when no external temperature 
gradient is present. The gas is Argon, the pressure 48 Pa, 
the RF amplitude 45V peak to peak, the discharge power 
17 mW. The number of particles is about 10^, the field of view 
32x43mm. The lowest 5 lattice planes of large interparticle 
distance consist of agglomerates. 

For the thermophoretic measurements the predeter- 
mined gradient of 1170 K/m was applied, the plasma 
was started and the discharge power was set to a medium 
value of 20 mW. Particles were injected until the plasma 
appeared to be saturated. To accomplish a stable dis- 
charge in the quasi 0-g case, higher RF amplitudes and 
powers are required. A survey was measured at three 
power levels (40, 57, 72 mW) and four pressures (14, 24, 
38, 46 Pa). 

Fig. m and ^ are taken out of the survey in order to 
visualize the stress-free complex plasma under quasi 0-g 
and in order to show a saturation of thermophoresis in 
the presence of a large particle number combined with a 
small working pressure. Fig. |^ shows how the particle 
cloud fills the whole discharge volume apart from space 
charge regions near to the electrodes and the central re- 
gion where a void is formed ^ . The reason for no lattice 
structure in Fig. ^ could be a residual convective motion 
in the order of Imm/s which is not visible in a single 
video frame. The void boundary appears well defined 
with a slight enhancement of particle density behind the 
boundary. Using A as given in table 1 we calculate a 
total heat flow from the lower to the upper electrode of 
24 mW. On the other hand the RF power dissipated in 
the plasma is roughly 57 mW. Most of this power will 



be dissipated by ionisation and heating of the electrodes 
by ion impact. Only a few mW will heat the neutral 
gas directly, hence we expect that the temperature pro- 
file applied from outside for thermophoresis remains rea- 
sonably homogeneous for the discharge parameters used, 
although a small amount of Joule heating in the center 
of the chamber might contribute to the void formation. 
If we compare our experimental condition to the one re- 
ported by Jellum et al. [Q, we find from the amplitudes 
and the power levels that they drive their discharge at 5 
to 10 times the power density that we use. If we increase 
the discharge power by 10, we find that the particles are 
driven out of the bulk plasma by both the ion drag force 
and intrinsic thermophoresis. 




FIG. 3: Side view of a complex plasma at quasi 0-g accom- 
plished by thermophoresis. The temperature gradient applied 
from outside is 1170 K/m. The peak to peak amplitude is 
82V, the discharge power 57 mW, the pressure 46 Pa and the 
number of particles injected in the order of 1 million. 

Obviously thermophoresis will not support an unrea- 
sonably large total mass of particles. A saturation in the 
thermophoretic force becomes visible first at low work- 
ing pressure. Fig. is a measurement under identical 
conditions as Fig. ^, however, at a working pressure 
of 14 Pa. The particles settle down towards the lower 
electrode. This effect can be reversed by increasing the 
pressure again to 46 Pa. Out of a number of possible 
mechanisms our data as well as rough estimates identify 
radiation cooling (Brattli and Havnes |]l6j) as the best 
explanation. By thermal radiation the particles will ap- 
proach everywhere a temperature roughly equal to the 
average temperature of the chamber surfaces. Hence for 
a large particle density and good thermal coupling be- 
tween the particles and the gas the temperature gradient 
in the dust region will be reduced and the gradient in the 
dust free region will become correspondingly larger. A re- 
duced temperature gradient in the dust region weakens 
the thermophoretic force. Gravity gains the overhand 
and the dust settles down towards the lower (heated) 
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electrode as observed. At larger pressure the decreasing 
mean free path of the neutral atoms reduces the ther- 
mal coupling between the gas and the particles, causing 
the reversal of this effect. A quantitative understanding 
of this effect requires detailed investigations beyond the 
scope of the present Letter. 




FIG. 4: Side view under identical conditions as Fig. ^ with 
only the gas pressure reduced from 46 Pa to 14 Pa. As the 
thermophoretic force becomes saturated particles settle down 
below the void. Interesting is a strong enhancement of particle 
density behind a very sharp void boundary. 

At ambient pressure (10'^ Pa) gravity driven convection 
would obstruct thermophoresis completely. At working 
pressures between 46 and 14 Pa convection is 3 orders of 
magnitude weaker. The driving force of vorticity depends 
on the temperature difference, the free fall acceleration 
and the mass density of the gas. The driving force hence 
decreases in proportion to pressure whereas the counter- 
acting friction on the chamber walls remains unchanged 
as the viscosity of gases is pressure independent. 

In conclusion, the thermophoretic levitation of mi- 
croparticles in a complex plasma, can give interesting 
insights into complex plasmas and be used as a promis- 
ing tool complementary to microgravity experiments. As 
an example we discussed radiation cooling by the parti- 
cles, leading to the observed particle density dependent 
reduction of the temperature gradient. Other important 
exampels are the investigation of the void and vortices 
formation, also observed under microgravity conditions. 
Furthermore we present an easy to use relation for the 
thermophoretic force together with the relevant gas ki- 
netic data. Technological applications of thermophoresis 



could be size selective particle production in a reactive 
plasma and surface plating of particles. In semiconductor 
production where nanometer sized particles are a byprod- 
uct in some processes, controlled thermophoresis may be 
used to keep particles away from the substrate where they 
can diminish the yield, as already suggested by Jellum et 
al. [0 , or to incorporate them in a controlled way where 
they are beneficial, e.g. for amorphous solar cells [ p7[ . 
Controlled thermophoresis could be interesting for "pick 
and place" in context with hybrid integration [ p^ . 
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